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Abstract. We report a detection in the interstellar medi- 
um of an infrared transition within the electronic ground 
state of the deuterated hydrogen molecule, HD. Through a 
deep integration with the Short- Wavelength-Spectrometer 
(SWS) on board the Infrared Space Observatory (ISO), 
the pure rotational v — — R(5) line at 19.43 /im was 
detected toward the Orion (OMC-1) outflow at its bright- 
est H2 emission region. Peak 1 . The ~ 20" beam-averaged 
observed flux of the hue is (1.84 ± 0.4) x 10~^ erg cm~^ 
s~^sr~^. Upper flux limits were derived for sixteen other 
rotational and ro-vibrational HD lines in the wavelength 
range 2.5 to 38 fim. 

We utilize the rich spectrum of H2 lines observed at 
the same position to correct for extinction, and to de- 
rive a total warm HD column density under the assump- 
tion that similar excitation conditions apply to H2 and 
HD. Because the observed HD level population is not 
thermalized at the densities prevailing in the emitting 
region, the total HD column density is sensitive to the 
assumed gas density, temperature, and dissociation frac- 
tion. Accounting for non-LTE HD level populations in a 
partially dissociated gas, our best estimate for the total 
warm HD column density is A^(HD) = (2.0 ± 0.75) x 
10^^ cm~^. The warm molecular hydrogen column density 
is (2.21 ±0.24) X 10^^ cm~^, so that the relative abundance 
is [HD]/[H2] = (9.0 ±3.5) x 10~^. 

The observed emission presumably arises in the warm 
layers of partially dissociative magnetic shocks, where HD 
can be depleted relative to H2 due to an asymmetry in the 
deuterium-hydrogen exchange reaction. This leads to an 
average HD depletion relative to H2 of about 40%. Cor- 
recting for this chemical depletion, we derive a deuterium 
abundance in the warm shocked gas, [D]/[H]= (7.6 ± 2.9) 
xlO"^. 



Send offprint requests to: F. Bertoldi 
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struments funded by ESA Member States (especially the PI 
countries: France, Germany, The Netherlands and the United 
Kingdom) and with the participation of ISAS and NASA. 



The derived deuterium abundance is not very sensitive 
to the dissociation fraction in the emitting region, since 
both the non-LTE and the chemical depletion corrections 
act in oppositite direction. Our implied deuterium abun- 
dance is low compared to previous determinations in the 
local interstellar medium, but it is consistent with two 
other recent observations toward Orion, suggesting that 
deuterium may be significantly depleted there. 

Key words: Shock waves - ISM: abundances - ISM: in- 
dividual objects: Orion Peak 1 ~ ISM: molecules - Cos- 
mology: observations - Infrared: ISM: lines and bands 



1. Introduction 

Deuterium is an important clue to the physics of the Big 
Bang. Its creation rate during primordial nucleosynthe- 
sis depended strongly on the number ratio of photons to 
baryons, a quantity not at all well-known but crucial for 
a correct description of the earliest events (e.g., Wilson & 
Rood 1994; Smith et al. 1993). Ever since, conditions have 
not been right to add further deuterium to the primor- 
dial production; neither nuclear fusion processes nor the 
spallation of heavier nuclei by energetic cosmic rays can 
augment the original abundance (although stellar flares 
were suggested by MuUan & Linsky [1999] to produce deu- 
terium). The deuterium abundance in fact decreases con- 
tinuously as deuterium is burned up in stars. The present 
day deuterium abundance in the interstellar medium pro- 
vides a lower limit to its primordial value, and it reflects 
the history of stellar reprocessing of the gas; measurements 
of its spatial variation may shed light on the star forma- 
tion history in a given region (e.g., Tosi 1998, Tosi et al. 
1998). 

Deuterium and D-bearing species were detected in gi- 
ant planets (Encrenaz et al. 1996; Feuchtgruber et al. 
1997), in cosmic rays (Beatty et al. 1985; Casuso & Beck- 
man 1997; Geiss & Gloeckner 1998), in the local ISM (Lin- 
sky et al. 1993, 1995; Linsky 1998; Piskunov et al. 1997; 
Helmich et al. 1996; Jacq et al. 1993; Turner 1990), and 
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in extragalactic sources (Songaila et al. 1994; Tytler et 
al. 1996; Buries & Tytler 1998). Through recent Di Lya 
observations of the local ISM, McCuUough (1992), Lin- 
sky et al. (1993, 1995), Piskunov et al. (1997), and Bring 
et al. (1997) derived a deuterium abundance [D]/[H] ~ 
1.5 X 10^^ along the line of sight toward nearby bright 
stars. These observations are hampered however by the 
possible confusion with the Hi Lya line, which is only 
~80 km apart. High spectral resolution UV absorp- 
tion observations of DI Lyman S and e were recently per- 
formed by Jenkins et al. (1999) and Sonneborn et al. (in 
prep.) along three lines of sight, which yield significant 
[D]/[H] variations between 7.4 x 10"^ and 2.1 x 10"^ Ob- 
servations of the hyperfine transitions of Di and Hi at 
radio frequencies (Chengalur et al. 1997), or of the H2 
and HD molecules, provide further means to sample the 
deuterium abundance. The [HD]/[H2] ratio is subject to 
variations due to small differences in chemical reaction 
rates in high temperature molecular gas as in shocks (see 
Sect. 3^) and photodissociation regions. To date, rota- 
tionally excited HD has been detected with ISO from giant 
planets (Feuchtgruber et al. 1997), and in the interstellar 
medium through UV-observations by Wright & Morton 
(1979). Accounting for the three lowest rotational levels 
of HD Wright & Morton found [HD]/[H2] < 6 x lO"'^ in 
the cold molecular gas toward ( Oph. D-bearing species 
such as HDO, CH2DOH, DCO+ or DCN (Jacq et al. in 
prep.), are observed at radio- frequencies in the ISM, but 
due to fractionation the deuterium abundance cannot be 
deduced accurately from these observations. 

The infrared emissivity of HD has been modeled and 
predicted for photodissociation regions (Sternberg 1990 - 
not including the reaction H-I-HD H2-I-D) and magnetic 
(C-type) shocks (Timmermann 1996), where latter calcu- 
lations include deuterated species in the chemical reaction 
network. But even state-of-the-art IR-spectrometers had 
been unsuccessful in detecting HD emission in the ISM, 
because of the low HD abundance, and because the pure 
rotational and rotation-vibrational lines appear at wave- 
lengths that are affected by strong telluric interference. 

ISO now for the first time opened the skies to a success- 
ful search for HD emission. Wright et al. (1999) detected 
the emission of HD 0-0 R(0) toward the Orion Bar pho- 
todissociation region, and derive [HD]/[H2]= (2.0 ±0.6) x 
10^^. A measure of the deuterium abundance in PDRs is 
complicated by the fact that the HD dissociation front is 
located deeper into the molecular cloud than that of H2 . 
The average excitation of HD may therefore be lower than 
that of H2 , making it difficult to derive column densities 
referring to the same regions. In regions of pure coUisional 
excitation such as shocks, this problem should not occur. 

We here report HD observations toward the bright- 
est H2 emission region in the sky, Peak 1 in the Orion 
molecular outflow that surrounds several deeply embed- 
ded far-infrared sources (Genzel & Stutzki 1989; Menten 
& Reid 1995; Blake 1997; van Dishoeck et al. 1998; Stolovy 



et al. 1998; Schultz et al. 1998). Our observations, which 
we discuss in Sect. 2, were part of a line survey of shocked 
molecular gas, and of a program to investigate the oxygen- 
chemistry in the warm molecular interstellar medium with 
the Short- Wavelength-Spectrometer (SWS). In Sect. 3 we 
derive an HD column density from the flux of one detected 
HD line, making use also of extinction and excitation mea- 
surements from a large number of H2 lines which we de- 
tected in related observations. We then estimate the deu- 
terium abundance in the shocked, warm gas of the Orion 
outflow. In Sect. 4 we summarize our results. 



2. Observations 

Orion Peak 1 was observed in the SWS 01 (full grating 
scan) and SWS 07 (Fabry-Perot) modes of the short wave- 
length spectrometer (de Graauw et al. 1996) on board 
ISO (Kessler et al. 1996) on October 3, 1997, and in the 
SWS 02 (« 0.01 A range grating scan) mode on Septem- 
ber 20, 1997 and February 15, 1998. Figure illustrates 
the various aperture orientations with respect to the H2 
1-0 S(l) emission observed with NICMOS on the HST 
(Schultz et al. 1998). The full 2.3 to 45 SWS 01 spec- 
trum was recorded in its slowest mode with the highest 
possible resolution. A preliminary reduction of this spec- 
trum was presented by Bertoldi (1997). Table |l| summa- 
rizes the HD line observations. The H2 lines will be dis- 
cussed in more detail in a forthcoming article (Rosenthal 
et al. in prep.). 

The data reduction was carried out using standard Off 
Line Processing (OLP) routines up to the Standard Pro- 
cessed Data (SPD) stage within the SWS Interactive Anal- 
ysis (lA) system. Between the SPD and Auto Analysis Re- 
sult (AAR) stages, a combination of standard OLP and in- 
house routines were used to extract the individual spectra. 
The in-house routines included an interactive dark-current 
subtraction for individual scans and detectors as well as 
for the removal of fringes. The flux calibration errors range 
from 5% at 2.4 /xm to 30% at 45 ^lm.. (SWS-Instrument 
Data User Manual, version 3.1). The statistical uncertain- 
ties derived from the line signal to noise ratio are for most 
detected lines smaller than the systematic errors due to 
flux calibration uncertainties. 



2.1. Line Fluxes 

The spectra at the expected wavelengths of seventeen HD 
transition lines are shown in Figs. ^ and |^. The wave- 
lengths (Table |l|) of the w = 0-0 pure rotational HD tran- 
sitions were adopted from Ulivi et al. (1991), while those 
for the rotation-vibrational lines were computed from the 
rotational (Essenwanger & Gush 1984) and vibrational 
constants (Herzberg 1950). Only one of the seventeen lines 
sought was detected, the 0-0 R(5) (i.e. J = 6 ^ 5) dipole 
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Fig. 1. The various SWS 
apertures of our ISO obser- 
vations of the Orion out- 
flow, overlaid on a NIC- 
MOS H2 1-0 S(l) map 
kindly provided by E. Erick- 
son and A. Schultz (Schultz 
et al. 1998). White dots 
and the central white patch 
are continuum ghost-images 
of stars and of Orion- 
BN, respectively The ISO 
aperture was centered on 
05'^32™46^27, -05° 24' 02" 
(1950). Its size corresponds 
to 14" X 20" for wavelengths 
smaller 12 /um, 14" x 27" at 
12 to 27.5 Aim, 20" x 27" at 
27.5 to 29 /im, and 20" x 33" 
at 29 to 45.2 ^m. 



transition at an expected wavelength of 19.4305 /irry (ob- 
served at 19.4290 /xm, see Fig. ||). Although the detection 
may seem marginal, the line does appear in two indepen- 
dent observations. It is apparent in the fringed data, and 
after defringing, it stands out well above the noise in the 
central wavelength range of the coadded SWS 02 scans. A 
beam-averaged flux of (1.84±0.4) x 10~^erg s~^cm~^sr~^ 
was derived from a continuum-subtracted integration over 
the feature; the error derives from the line's S/N= 4.5 (the 
RMS noise being evaluated within ±500 km s~^) plus an 
estimated 11% flux calibration uncertainty. The line width 
(FWHM) is 134 km s~^, which agrees with the instrumen- 
tally expected width for extended objects, 130 km s^-'^, and 
is not sufficient to resolve the emission, which should have 
a velocity dispersion similar to that observed in II2 1-0 
S(l), which has a FWHM of « 50 km s"\ with emis- 
sion from -100 to -1-100 km s~^ (Chrysostomou et al. 
1997; Stolovy et al. 1998). The line center is positioned 
at i^hoiio — —23 km s~^, which is within the range of line 
center velocities, Whoiio — —38 to -1-41 km s~^, observed 
in the H2 1-0 S(l) emission (Chrysostomou et al. 1997). 
The high spectral resolution and sensitivity of the ISO 
SWS around 19 /j,m made — R(5) the most promising 
line for a detection of HD in Peak 1. The stronger lines 
from lower rotational states suffer from the rapidly rising 



^ The experimental value for the R(5) transition wave- 
length is 19.4305 ± O.OOOlpim, whilst the calculated value 
is 19.431002 ± O.OOOOOS^tm (Ulivi et al. 1991 and references 
therein). The discrepancy is yet unexplained. 



continuum level at longer wavelengths and the resulting 
strong fringing. 

The available integration time was insufficient to de- 
tect two other HD lines which we tried to observe in the 
SWS 02 mode. We derived upper flux limits for a total of 
sixteen HD lines from the noise level near the expected line 
wavelengths (Table |^). The peak-to-peak noise envelope 
approximately corresponds to a 3ct flux density disper- 
sion, and we assumed that a line with a peak flux density 
of 3a would have stood out clearly enough to be detected. 
As upper flux limits we therefore adopted the 3f7 flux den- 
sity noise level times the FWHM of neighboring H2 lines, 
except for 0-0 R(7), R(8), R(9), for which we adopted the 
linewidths expected from instrumental resolution. 



3. Extinction correction and H2 excitation 

The column density of molecules in a particular rotation- 
vibrational level {v, J) is computed from the observed line 
flux, Iobs{v,J v',J'), of a transition to a lower state 
(w', J'), times the wavelength, A, divided by the radiative 
transition rate (Einstein) coefficient^ for this transition, 



The Einstein coefficients for H2 and HD are taken from 
Turner et al. (1977), confirmed by Wolniewicz et al. (1998), 
and Abgrall et al. (1982), respectively. The transition energies 
for H2 were computed from the level energies that were kindly 
provided by RouefT (1992, private communication). 
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Table 1. Summary of ISO-SWS observations of HD lines at Orion Peak 1. 



transition 


A 


b 

integration peak 


Air ^ 


A 


P -r 


-^obs 




A\ 


AT 






sec Jy 


km S 


s 




±v 


erg s '^cm ^sr ^ 


mag 


cm 


0-0 R(2) 


37.7015 


< 192 


300 


1.72 




765.9 


<9.9 


-4) 


0.09 


< 1.5(17) 


0-0 R(3) 


28.5020 


< 100 


350 


4.11( 


-6) 


1270.7 


< 7.9 


-4) 


0.30 


< 4.6(16) 


0-0 R(4) 


23.0338 


< 35 


350 


7.91( 


-6) 
'-'1 


1895.3 


< 6.0 


-4) 


0.51 


< 1.8(16) 


0-0 R(5) 


19.4305 


3010 2.37 


134 


1.33( 


-5) 


2635.8 


(1.84 ± 0.4) 


-5) 


0.60 


(3.0± 1.1)(14) 


0-0 R(6) 


16.8940 


< 8.5 


230 


2.03( 


-51 


3487.5 


< 1.3 


-4) 


0.53 


< 1.1(15) 


0-0 R(7) 


15.2510 


1110 < 5.0 


185 


2.88( 


-5) 


4445.3 


< 6.8 


-5) 


0.41 


< 3.3(14) 


0-0 R(8) 


13.5927 


< 10 


215 


3.87( 


-5) 


5503.8 


< 1.8 


-4) 


0.37 


< 5.6(14) 


0-0 R(9) 


12.4718 


< 8 


230 


4.97( 


-5) 


6657.5 


< 1.7 


-4) 


0.54 


< 4.3(14) 


1-0 P(4) 


3.0690 


< 0.2 


280 


7.37( 


-6) 


5958.5 


< 2.8 


-5) 


1.11 


< 2.0(14) 


1-0 P(3) 


2.9800 


< 0.1 


270 


1.09( 


-5) 


5593.5 


< 1.4 


-5) 


1.08 


< 6.4(13) 


1-0 P(2) 


2.8982 


< 0.25 


300 


1.64( 


-5) 


5348.7 


< 3.9 


-5) 


0.94 


< 1.0(14) 


1-0 P(l) 


2.8225 


< 0.4 


305 


3.23( 


-5) 


5225.7 


< 6.6 


-5) 


0.80 


< 7.6(13) 


1-0 R(0) 


2.6900 


< 0.2 


330 


1.72( 


-5) 


5348.7 


< 3.7 


-5) 


0.70 


< 7.0(13) 


1-0 R(l) 


2.6326 


< 0.3 


330 


2.51( 


-5) 


5593.5 


< 5.7 


-5) 


0.70 


< 7.3(13) 


1-0 R(2) 


2.5811 


< 0.4 


245 


3.22( 


-5) 


5958.5 


< 5.8 


-5) 


0.72 


< 5.7(13) 


1-0 R(3) 


2.5350 


1200 < 0.2 


215 


3.91( 


-5) 


6441.1 


< 2.6 


-5) 


0.74 


< 2.1(13) 


1-0 R(4) 


2.4943 


< 0.35 


235 


4.60( 


-5) 


7037.7 


< 5.0 


-5) 


0.76 


< 3.5(13) 



° expected wavelengths. 0-0 R(5) is observed at 19.4290 /xm. 

' on target integration time for SWS 02 observations; the SWS 01 2.4 - 40 fim scan took 6538 sec. 

FWHM of observed HD line, or of neighboring H2 lines. 

upper level column, corrected for extinction. 
" numbers in brackets denote powers of ten. 



A{v,J — > w', J'), times a correction for extinction along 
the line of sight: 



^ ^ 47r A Iohsiv,J-^v',J') 
"■^ he A{v,J^v',J') 



10' 



(1) 



where A\ is the effective extinction at wavelength A. All 
H2 and HD transitions are optically thin. Previous mea- 
surements of the Peak 1 molecular emission estimated K 
band (2.12 jj,m) extinctions between 0.5 and 1 magnitude 
(Everette et al. 1995). To apply an extinction correction to 
the observed H2 and HD lines requires knowledge of the 
"extinction law" A\/Ak between 2.4 and 40 /xm, which 
especially above 5 ^m is observationally not well con- 
strained (Draine 1989). The near-IR extinction is usually 
approximated to follow a power law (see Fig. ^), 



Ax^Ak ( a / 2.12 ^m 



for A < 6 /im, 



(2) 



with an increase in opacity beyond 6 /j,m due to stretch 
and bend mode resonances in silicate grains (Draine & 
Lee 1984), the shape and depth of which are yet poorly 
understood. A recent ISO study of the extinction toward 
the W51 Hll region, e.g., finds Aig/Aio - 0.52 ±0.1 and 
Aiq/Ak ^ 0.57 ± 0.1 (Bertoldi et al. in prep.). Draine & 
Lee (1984) suggested Ais/Ag,^ ~ 0.40, and observations of 
circumstellar dust emissivities suggest values for this ratio 
between 0.35 (Pegourie & Papoular 1985) and 0.5 (Volk 
& Kwok 1988). 



The effective extinction toward the emitting gas in 
Peak 1 may not necessarily follow an average interstel- 
lar extinction law, since the emitting and absorbing gas 
may be mixed. We therefore tried to estimate the effec- 
tive extinction as a function of wavelength toward Peak 
1 from the observed emission of Peak 1 itself. The differ- 
ential extinction between two wavelengths can be derived 
from a comparison of H2 line fluxes of transitions arising 
either from the same upper level, or from neighboring ther- 
malized levels with a well determined relative excitation. 
Since the warm HD and H2 are likely to be well mixed, 
we will use the extinction toward the H2 to deredden also 
the HD line intensities. 



The extinction-corrected H2 column density distribu- 
tion serves as a "thermometer" that probes the excitation 
conditions in the emitting gas. Since we detected only one 
HD line, we cannot determine the HD excitation directly 
-this would require at least two lines of transitions arising 
from different upper levels. Instead, we rely on the reason- 
able assumption that HD is subject to the same excitation 
conditions as H2. We can thereby derive both molecules' 
total column densities in the warm, emitting gas, and from 
that, their abundance ratio. 
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Fig. 2. Spectra of HD non-detections. Three of the spectra were obtained in the SWS 02 mode with ISO, all others 
were taken in the SWS 01 mode. Neighboring H2 and HI lines are identified at their expected positions. 



3.1. 2.4 - 6 fim extinction 

It is useful to plot the observed H2 or HD column den- 
sities, Nyj, divided by the state's statistical weight]^ gj, 
against the state's energy, Eyj. Examining this "excita- 
tion diagram," which we first plotted from the not yet 
extinction-corrected H2 line intensities, we found no indi- 
cation up to Eyj/k « 40,000 K of fluorescent excitation, 
or for deviations from the statistical ortho- to para-H2 
abundance ratio of 3 (details will be given by Rosenthal 



^ The statistical weights are gj = 2J+1 for para-H2 (even J) 
and gj = 3(2J-f 1) for ortho-H2 (odd J). HD has non-identical 
nuclei (and a small dipole moment ~ 10"'* Debye). Therefore, 
unlike H2, which exists in the form para-H2 and ortho-H2, such 
a difference does not exist for HD, and gj = 2J + 1. 



et al. in prep.). Instead, the populations appear both in 
thermodynamic and statistical equilibrium, which means, 
they depend only on the level energy and degeneracy, not 
on the vibrational quantum number. As a consequence, 
Nyj/gj is a smooth function of the level energy. 

In contrast, H2 that is fluorescently excited generally 
displays level populations that are not in vibrational and 
rotational LTE, and also show ortho-to-para column den- 
sity ratios between neighboring states that are smaller 
than 3 (see e.g., Draine & Bertoldi 1996, Bertoldi 1997). 
Timmermann (1998) predicted that in low-velocity shocks 
the ortho-to-para ratio can also be significantly lower, 
which was confirmed in recent ISO-SWS observations by 
Neufeld et al. (1998), who found a value of ~ 1.2 in the 
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Fig. 3. Befringed AOT 02 spectra of two 0-0 R(5) obser- 
vations performed at the labeled dates, and the merged 
spectrum. The line is visible even in the original fringed 
data of both observations. Continuum levels differed in 
both datasets (due to differing orientations of the aper- 
ture with respect to Orion-BN), so that a constant was 
added to one before averaging them, and in order to dis- 
play them together. The rising noise near the scan edges 
is expected and due to the lower detector coverage. In the 
merged spectrum, the line stands out with a peak flux den- 
sity to RMS noise ratio of 4.5. Whciio refers to the expected 
line wavelength of 19.4305 /im. 

outflow HH54. Peak 1 however shows no deviations from 
statistical equilibrium. 

We were able to obtain the most reliable line fluxes for 
the H2 w = — and ?; = 1 — lines between 2.4 and 6 /im 
that have upper level energies E^j jk = 5000 — 16, 000 K. 
Excluding lines in the water ice absorption feature be- 
tween 2.8 and 3.3 ^m, we selected thirty-two high S/N 
lines in this wavelength and energy range to constrain the 
power law part of the extinction curve. Assuming that 
the H2 columns N^j/gj vary smoothly with E^j, we fit 
a second order polynomial to the observed column den- 
sities, corrected for extinction following Eq. (1) with 
and e as free parameters. We searched for values of Ayl 
and e which minimize the dispersion of N^^j/gj. Within 
the range e = 1 — 2, Ayi is well constrained to 1.0 ± 0.1 
magnitudes. It was not possible to narrow down the value 
of e further, so we adopted e — 1.7, which was suggested 
by previous observational studies (see Draine 1989; Brand 
et al. 1988 used — 0.8 and e — 1.5). The second order 
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Fig. 4. Adopted extinction curve (solid line), includ- 
ing water ice and silicate features. Dashed lines show 
curves with different relative silicate features strengths, 
Aig/AiQ = 0.35 and 0.52. The dash-dotted curve assumes 
narrower 9.7 fxm feature. Wavelenghts of those H2 lines 
used to constrain the extinction curve are marked - but 
these are not datapoints! 

fit to N^j/gj versus E^j/k = 1000 T3 K for Ak = 1.0 and 
e = 1.7 is (see Fig. 1) 

\og{N^j/gj) = 18.88 - 0.402 Ta-f 0.0092 T^, (3) 

for 5800 K < E^j/k < 17,000 K, and has a fit quality 
— 0.13; the uncorrected {Ak — 0) dispersion is = 
0.39. 

We should note that our extinction curve between 4 
and 7 fim is not constrained well enough to address the 
claim by Lutz et al. (1997) that the extinction curve flat- 
tens and lacks the presumed minimum at 7 /im. 

3.2. Mid-IR extinction 

Having constrained the extinction curve below 6 /xm and 
the H2 excitation from 5800 - 17,000 K, wc now estimate 
the strength of the 10 /j,m silicate feature with the ob- 
served H2 w = 0, J < 8 column densities. The intensity 
of the 0-0 S(3) line at 9.66 /im gives a column density 
of the w = 0, J = 5 level that is significantly below that 
expected from an interpolation between the J — i (0-0 
S(l) 17.0 /im) and J = 8 (0-0 S(6) 6.1 /im) levels. The 
observed J = 3 to J = 7 level populations are all af- 
fected by silicate absorption. We estimate the shape of 
the 9.7 /im and 18 /im silicate features from the calcula- 
tions of Draine & Lee (1984), and adjust the strengths of 
the features individually (Fig. Assuming A^i = 1.0 mag 
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Fig. 5. H2 excitation diagram: the dots, diamonds, and 
squares denote dereddened and beam-averaged column 
densities of transitions in the v = 0-0, 1-0, and 2-1 bands, 
respectively. The solid line represent the least squares fit, 
Eqs. (||) and (Q), which apply for E.^j/k larger and smaller 
than 5800 K, respectively. 



and tuning ^9.7 such that the J — b level column density 
smoothly follows that of the J = 3 — 8 levels (Fig. |5| shows 
the resulting dereddened columns) results in a possible 
range A9.7 — (1.3 — 1.5) mag, if we allow for a variation in 
Aig/Aio — 0.35 — 0.52, and a 15% uncertainty of the J — 5 
column due to the difference between two measurements 
of the 0-0 S(3) line in different AOTs. When we narrow 
the 10 /im feature below reasonable width estimates, the 
extinction correction for the J — A and J = 6 levels de- 
creases significantly, thereby lowering the dereddened flux 
in the 0-0 S(3) line and decreasing the required extinction 
to Agj « (1.1 — 1.3) mag. 

Considering all uncertainties we estimate ^9.7 ~ 
(1.35 ± 0.15) mag, and with Aig/Aw « 0.35 - 0.52, we 
find Aig Ri (0.61 ± 0.15) mag. Our value for Aiq/Ak = 



(0.61 ± 0.15) is in good agreement with the 0.57 ± 0.1 
found toward W51. 

The extinction correction for the detected HD line we 
then estimate as lO^-^^^s ~ 1.75 ± 0.25. The dereddened 
column densities of the H2 w = 0, J < 8 levels are fit by 



log {Nyj/gj 



20.17- 0.765 Tg + 0.0344 To^ 



(4) 



At upper level energies Ej^/k = 0, 2636 K (HD J = 6) 
and 5800 K this corresponds to excitation temperatures 
Te^ ~ 570, 740 and 1190 K, respectively. 



3.3. Water ice feature 

The rotation- vibrational HD transitions 1-0 P(2), P(3), 
and P(4) fall into the water ice absorption feature at 
3.05 ±0.25 fim. To correct these line flux limits for extinc- 
tion we estimated the depth and width of the feature from 
the apparently enhanced (over the power law, Eq. [^) ex- 
tinction of five H2 lines with 6000 K < Eu/k < 16, 000 K 
and 2.8 /im < A < 3.3 /im. The Gaussian 

-[(A-3.05/im)/(%/2-0.15Aim)] 



AA,,,{X) w 0.58 e- 



mag 



(5) 



approximately fits the additional extinction noticeable for 
the five H2 lines, and was used to correct the rovibrational 
HD line fluxes. 

3.4. Total column 0/H2 andUD 

How can we derive the total HD column density with only 
one measured high-excitation level? We can make use of 
the observed H2 excitation, and estimate that the excita- 
tion conditions - i.e., the fractions of the total gas col- 
umn density that are at a particular temperature and 
density - are the same for HD and for H2. We assume 
that at least up to the excitation energy of HD J — 6 
{Eoe/k = 2636 K), the H2 level populations are thermal- 
ized, and thereby reflect the kinetic temperature distribu- 
tion of the gas. The assumption of thermalized level pop- 
ulations is supported by the lack of deviations from vibra- 
tional degeneracy or non-statistical ortho-to-para ratios, 
and by detailed non-LTE calculations of the H2 level popu- 
lations, which show that deviations from LTE are small at 
J <8{E/k< 5800 K) for densities higher than 10^ cm'^ 
and temperatures above 600 K (Draine & Bertoldi, un- 
published) . 

Because they are permitted dipole transitions, the ra- 
diative decay rates of HD are much higher than those of 
H2 levels at comparable energy. This results in "critical" 
densities (i.e. densities above which a level is thermalized 
for a given temperature) which are higher for HD than 
for H2. If the gas density was high enough to thermalize 
the HD levels, then HD would show the same level exci- 
tation as H2 : normalizing the H2 populations (|^) with the 
measured HD J — 6 level, we would then expect the HD 
populations to follow 



15.14 - 0.765 T3 + 0.0344 T| 



(6) 
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Fig. 6. HD excitation diagram. Pure rotational transi- 
tions are denoted by dots, and v = 1-0 transitions by 
diamonds. The line represents the fit Eq. (^. The error of 
the 0-0 R(5) line is computed from spectral noise (22%) 
and uncertainties in the flux calibration (11%) and the 
extinction at 19.4 /im (~14%). 

for T < 6000 K. The measured upper limits to sixteen HD 
level column densities are consistent with this distribution 
(Fig. I). 

However, deviations from LTE are important at the 
expected density of 10^ — 10^ cm~'^ and temperatures of 
600-1000 K in the shocked gas of the OMC-1 outflow. Let 
n„,//nLTE,ti,/ be the actual non-LTE population in level 
{v, J), divided by its LTE value. The total warm HD col- 
umn density is given by the sum over all level column 
densities 



7V(HD) 



»^ LTE, 06 
no6 



E 

vj 



9J 



(2J+1) 



(7) 



where Nyj/gj is given by Eq. (^). 

To compute the deviations from LTE of the HD level 
population we solved the equations of statistical equilib- 
rium (see Timmermann 1996), including radiative decay 
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Fig. 7. Deviations from LTE HD level populations plot- 
ted against the rotation quantum number J for different 
densities and kinetic gas temperatures, tiqj and nLTE,oj 
refer to densities in an individual v — state under non- 
LTE and LTE conditions. 



and excitation/deexcitation through collisions with He, 
H2, and H. The abundance of H varies throughout the 
shock, and we adopted [H]/[H2]= 1, which is typical for 
the dissociation fraction in the hot layers of a partially 
dissociative shock. In Fig. ^ we compare the resulting 
level population to that in LTE for gas volume densities 
of 3 X 10^ and 3 x lO^cm^^, and three different kinetic 
temperatures that should span the range typical of the 
observed, shocked gas. 

HD collisional (de)excitation rate coefficients with H2 
were computed properly only for pure rotational transi- 
tions up to J = 4 and temperatures up to 600 K and 
300 K, respectively (Schafer 1990). The rate coefficients 
for higher level transitions and temperatures were extrap- 
olated from the lower ones. The H-HD and H2-HD colli- 
sion rate coefficients were recently computed by Roueff & 
Flower (1999) and Roueff & Zeippen (1999), respectively. 

We find that for H2 densities below 10^ cm"'^ the HD 
level populations deviate significantly from LTE. At a den- 
sity of 3 X 10^ cm^^ and a temperature of 600 K, the pop- 
ulation of the J = 6 level is about a factor two below its 
LTE value. 

To assess the non-LTE effects on the derived total HD 
column density, we evaluated A^(HD) from Eq. (^ and 
plot this in Fig. || as a function of the gas density for 
three different temperatures. Only at H2 densities above 
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Fig. 8. Total HD column density computed from non-LTE 
level populations as a function of H2 density for three dif- 
ferent gas temperatures. The H2 density in the preshock 
gas has been estimated at (1 — 3.5) x lO^cm"'^, and the 
shocked gas temperatures at w 600 — 1000 K, so we esti- 
mate iV(HD) « (2.0 ± 0.75) X lO^^cm-^ 

10^ cm~'^ do the populations assume LTE, and the warm 
HD column is 

7V(HD)lte = (1.36 ± 0.38) X lO^'^cm-^ . (8) 

The error derives from the uncertainty of the J ~ 6 col- 
umn measured by the 0-0 R(5) line, and is a combination 
of spectral noise (22%), flux calibration uncertainty (11%), 
and the uncertainty of the extinction correction (14%). At 
lower densities, the total column density much depends on 
the gas temperature, and since H is the strongest collision 
partner of HD, on the dissociation fraction. 

The emission of the OMC-1 outflow had been mod- 
eled previously with C-shocks that propagate at veloci- 
ties of order 35 — 40 km s~^ into gas with densities of 
n(H2) = (1 - 3.5) X 10^ cni-3 (Draine & Roberge 1982; 
Chernoff et al. 1982; Kaufman & Neufeld 1996). The gas 
in such shocks reaches temperatures of order 1000 K, and 
remains at approximately constant density through the 
region where most of the H2 emission occurs. The low- 
est H2 levels we observed show excitation temperatures 
of 600-700 K, which probably reflects the kinetic temper- 
ature of much of the warm, emitting gas. Taking these 
temperatures and densities as estimates for the prevailing 
excitation conditions, we estimate that the total observed 
HD column density must be in the range 

iV(HD) = (2.0 ± 0.75) X 10^*^ cm"2 ^ (9) 

If we had neglected H-HD collisions, the HD level would 
be less thermalized, and our estimate for the total HD col- 



umn would rise to (3.5 ± 1.4) x 10^^ cm~^. However, we 
believe that much of the emission arises from partially dis- 
sociative shocks, so that collisions with neutral hydrogen 
should be important. 

We now compare the HD column to that of warm H2 . 
Summing over the H2 level populations given by the least 
squares fits, Eqs. (||) and (^, we find 

Ar(H2) = (2.21 ±0.24) x 10^^ cm'^ , (10) 

where the error reflects a maximum flux calibration uncer- 
tainty of 11% in the 7 to 19.5 /im range. By summing from 
J = 0, we extrapolated the observed H2 level populations, 
J > 3, to the unobserved levels J = — 2. We thereby ac- 
count only for the warm H2, not for the total H2 column 
along the line of sight, which includes over 10^^ cm~^ of 
cold gas in the molecular cloud that embeds the outflow. 
In this cold gas, which we are not concerned with, most 
H2 is in its ground states J = and J = 1, and does not 
affect our analysis of the warm outflow gas seen in the 
higher levels. 

Dividing the HD and H2 column densities, we derive a 
first estimate of the abundance ratio 

[HD]/[H2] = (9.0±3.5) X 10"^ (11) 

3. 5. Chemical depletion of HD in shocks 

The OMC-1 outflow emission arises from warm molecu- 
lar gas that is shock-heated by high-velocity ejecta which 
originates from one of the deeply embedded protostars 
in the vicinity of Orion IRc2. The emission is a mixture 
from fast, dissociative J-type shocks, in which the molec- 
ular emission comes from where molecules reform, and 
from slower (< 50 km s~^), partially- or non-dissociative 
C-shocks, in which the molecules radiate in a magnetic 
precursor where the temperature rises to its peak value 
(Fig. In J-shocks that propagate into a medium with 
density of order 10^ — 10® cm^'^, the fraction of the bulk 
motion energy radiated away in H2 lines is very small, of 
order 0.1% (Neufeld & Dalgarno 1989). In C-type shocks, 
however, about half the kinetic energy is converted to H2 
emission. For J-type shocks to dominate the H2 and HD 
emission, about 1000 times more energy would have to be 
dissipated in fast J-shocks than in the slower C-shocks. 
For an even distribution of magnetic field strengths and 
shock velocities up to ~ 100 km s~^, C-shocks therefore 
vastly dominate the molecular emission. The outflow emis- 
sion may therefore be well modeled as arising from such 
magnetic, partially-dissociative shocks. 

Most of the deuterium is locked in HD in the dense 
pre-shock molecular gas. In a high- velocity C-shock, both 
HD and H2 are either dissociated through collisions with 
ions that stream through the neutral gas at a speed of or- 
der the shock velocity, or through collisions with warm H2 
and H at temperatures in excess of 2500 K. In a partly dis- 
sociative C-shock, HD is depleted more than H2 because 
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with the code of Timmermann (1996, 1998). The pre- 
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7.5 cm~'^, the shock speed Us = 35 km s~-^, and Bq = 
700/iG. The depletion of HD relative to H2 reaches a max- 
imum value of 2.4 behind the hottest layer of the shock. 



atomic hydrogen can efficiently destroy HD through the 
reaction 

D + H2 t^HD + H + Ai/o, (12) 

where AHo/k = 418 K is the enthalpy difference. Figure H 
illustrates the density profiles for H2 and HD across such a 
C-shock. In the warmest layers, HD is depleted by up to a 
factor 2.4 relative to H2. In the cooling region, the atomic 
deuterium reacts stronger again with H2 to form HD, and 
the equilibrium of Eq. ( [l^ ) is shifted back towards HD. 
Averaged over the region where the temperature exceeds 
400 K, from where most of the observable emission arises, 
HD is depleted relative to H2 by a factor 1.67 in the par- 
ticular case we display in Fig. 

In C-shocks propagating at velocities below 25 km s^^, 
H2 is not dissociated significantly, and because of the low 
abundance of atomic hydrogen, the chemical depletion of 
HD is negligible. 



3.6. Deuterium abundance 

To derive the deuterium abundance in the warm, shocked 
gas we adopt the non-LTE level distribution of HD we 
computed above, and further account for the possibility of 
chemical depletion. With the column density ratio derived 
for temperatures of 600-900 K and densities of (1 — 3.5) x 
10^ cm^'^, we found 

[D]/[H] = 0.5[HD]/[H2] = (4.5 ± 1.7) x 10"^ , (13) 

and accounting for chemical depletion of HD by 1.67, the 
abundance is raised to 

[D]/[H] = (7.6 ± 2.9) X 10~^ , (14) 

where no error was added for uncertainties in the chemical 
depletion factor. 

The main uncertainty of our value appears to arise 
from the indirect measure of the HD excitation, and of 
the abundance of atomic hydrogen in the warm shocked 
region. Less dissociative shocks would require larger cor- 
rections for non-LTE level populations, which would raise 
the implied HD abundance: neglecting all H-HD collisions, 
we earlier derived iV(HD) « (3.5 ±1.4) x lO^^cm"^, which 
yields [D]/[H]= 7.9 x 10~^. Now chemical depletion should 
not occur in shocks with a low abundance of neutral hy- 
drogen, so that this values reflects the actual deuterium 
abundance. Interestingly, we find that the effects of non- 
LTE and of chemical depletion nearly cancel, so that the 
derived deuterium abundance turns out to be not very 
sensitive on how dissociative the shock is. 

4. Summary 

ISO for the first time enabled the detection in the inter- 
stellar medium of an infrared transition in the electronic 
ground state of deuterated hydrogen, HD. We here re- 
port the discovery of the w=0-0 R(5) line at 19.4290 fim 
with the ISO Short Wavelength Spectrometer in the warm, 
shocked molecular gas of the Orion OMC-1 outflow, at the 
bright emission "Peak 1." Upper flux limits for sixteen 
other HD lines were measured, all of which appear consis- 
tent with expectations when considering the observed 0-0 
R(5) Hue flux. 

A large number of H2 lines were detected (Rosenthal et 
al. in prep.) and utilized to analyze the HD observations. 
The near- and mid-infrared extinction toward the emit- 
ting region was derived by minimizing the dispersion in 
the observed H2 level column densities with respect to an 
LTE excitation model. Thereby we derive a near-infrared 
(K band) extinction of (1.0 ± 0.1) magnitudes, a 9.7 jim 
extinction of (1.35 ± 0.15) mag, and from an estimated 
range of A19/A9.7 = 0.35 — 0.52 we correct the HD 0-0 
R(5) flux for extinction by (0.61 ± 0.15) mag, i.e. a factor 
1.75 ±0.25. 

The dereddened H2 level populations served to esti- 
mate the excitation conditions in the gas. While H2 was 
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assumed to have thcrmalized level populations, those for 
HD were computed in detail by making use of the H2 ex- 
citation. Due to non-LTE effects at J > 3, the total warm 
HD column density was found to be sensitive to the gas 
density and temperature at the densities estimated to pre- 
vail in the shocked gas, n(H2) w (1 - 3.5) x 10^ cm~^. 

Our estimate for the observed warm HD column den- 
sity is iV(HD) = (2.0±0.75) x 10^'^ cm'^, and for the warm 
molecular hydrogen, iV(H2) = (2.21 ± 0.24) x 10^^ cm'^. 
Their relative abundance is therefore [HD]/[H2] = (9.0 ± 
3.5) X 10-6. 

We note that in high-velocity C-shocks, HD may be 

depleted relative to H2 because of an asymmetry (due 
to a small binding energy difference) in the deuterium- 
hydrogen exchange reaction HD + H t5 D + H2 . Estimat- 
ing that this lowers the warm HD column by about 40%, 
we derive a deuterium abundance in the warm shocked 
gas, [D]/[H]= (7.6 ± 2.9) x 10"^ 

If the emitting shocks were on average less dissocia- 
tive than assumed, the chemical depletion would be less 
pronounced. But at the same time, the lower H-HD colli- 
sion rate would enhance the HD J = 6 level population's 
deviation from LTE, to the effect that our implied total 
column of HD would increase. We estimate that the two 
effects approximately cancel, and that [D]/[H] remains at 
« 8 X 10^6, independent of how dissociative the shock 
actually is. 

The major uncertainties in our estimate for the deu- 
terium abundance arise from our indirect measure of the 
HD excitation, and since we must therefore rely on non- 
LTE excitation models, on the uncertainty of the preshock 
density and the abundance of neutral hydrogen in the 
shock. Future ground-based near-IR observations of ro- 
vibrational transition lines of HD could constrain the ef- 
fects of non-LTE and thereby narrow the error margins 
of the deuterium abundance. Detailed shock models of 
the recently resolved, multiple bow-shaped emission in the 
Orion outflow (Schultz et al. 1998) might also yield better 
estimates for the pre-shock densities and shock velocities, 
and thereby of the H2 dissociation and of the chemical 
depletion fraction of HD in such shocks. 

The deuterium abundance we find, (7.6±2.9) x 10"^, is 
lower than that derived through most DI absorption mea- 
surements in the local ISM, but it is consistent with that 
found recently by Jenkins et al. (1999) toward A Orionis, 
(7.4l};3) X 10"*^, and by Wright et al. (1999) toward the 
Orion Bar, (1.0 ±0.3) x 10~^. This could indicate that to- 
ward Orion, the deuterium abundance is indeed somewhat 
lower than on average. 
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